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ABS ted 


A system for PERT analysis using Monte Carlo methods to estimate 


project completion times was designed u tng a 2250 Display Unit 


operating on an IBM 360/67 computer. The user enters the PERT network 


and the distributions of activity times by using the light pen on the 
2250 Display Unit. 
After the computer performs the analysis, the probability of being 


on the critical path is shown for each arc in the network and a histo- 


gram of project completion times is displayed. The user may then 


modify the network and run another analysis. 

A discussion of the difference between this Monte Carlo PERT 
analysis and the usual PERT analysis and details of the implementation 
A user 


is included, along with operating instructions for the program. 


with no programming experience on the 2250 Display Unit can utilize the 


program by following the operating instructions. 
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IL. INTRODUCTION 





PERT (Program Evaluation and Review Technique) is a technique for 
determining the expected time of completion and the variance associated 
with this expected time for a project involving a number of activities 
which are interrelated. The usual PERT analysis, described in ref. l 
and 2, involves analyzing the network once using expected activity 
durations and computing an expected completion time, a variance on this 
completion time, and a critical path. 

A Monte Carlo analysis, discussed in ref. 3, involves analyzing the 
network a great many times using random numbers to draw activity duration 
times from the distribution of each activity's duration times. This 
method also computes an expected project completion time and a variance 
on this completion time. However, inst.ad of a single critical path, a 
criticality index is calculated for each activity. This is the ratio 
of the number of times an.arc was on the critical path to the total 


number of analyses. 


A. STATEMENT OF THE PROBLEM 

To analyze even a fairly simple PERT network by hand can be very 
time consuming. To analyze a network many times as is done in a Monte 
Carlo analysis would be : Herculean task if done by hand. At the same 
time, a student who is to be expected to use PERT networks in the 
future should know how a network is analyzed and have a good underst: :ding 
Be the meaning Of critical path, criticality indexes, and th  cffects 
on these from cha. zes to the network and the activity duration times. 


One way to accomplish these objectives is to combine the rapid 








processing of the computer with a visual display which allows the student 
to experiment and change the network. 

This paper describes a computer program, PERT-O-GRAPH, which allows 
a user to draw a PERT network by using a light pen on a graphic display 
unit. Distributions of activity duration times can be entered by 
drawing the distributions with the light pen or reading them from 
punched cards. A single analysis, called a "Mean Time Calculation” in 
this program, can be made which is similar to the usual PERT analysis. 
The user can then run an analysis using Monte Carlo technigues to 
estimate the distribution of project completion times. The program 
illustrates the changing of the critical path as activity durations 
vary, and the criticality indexes of each arc on the network after the 
analysis is complete. 

PERT-O-GRAPH has i cen designed to operate on an IBM 360/67 computer 
with an IBM 2250 Display Unit Model 1. After reading this paper a user 
should be able to enter and eadyee PERT networks even if he has no 
previous experience with the IBM 2250. No knowledge of computer pro- 


gramming techniques is necessary to use the program. 


B. ORGANIZATION OF THE PAPER 

Section II of this report describes the usual PERT analysis. A 
method for analyzing a PERT network and an example analysis are given. 
Section III explains the difference between the usual PERT analysis and 
a Monte Carlo anzlysis. Section IV discusses the intended use: of the 


program and describes the graphic displays and the hard copy output. 
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Section V contains operating instructions for the use of PERT-O-GRAPH. 
The use of the light pen, the programed function keys, and the alpha- 
numeric keyboard associated with the 2250 Display unit are described. 

section VI is an enol problem. Jt contains photographs of the 
graphic displays and illustrates how the user can change the PERT 
network. 

Section VII is a description of the program. Concluding remarks 


appear in Section VIII, followed by the program listing. 


at 








De Rea AY os 


In the literature dealing with PERT networks the following terms 
are commonly used, and these definitions will be used throughout this 
discussion. The arcs of a PERT network are called activities and have 
associated with them so:.e duration called ee The nodes of a network 
are called events, and they represent the completion of activities. 

The PERT network is a description of the relationships among the 
activities. These activities must be finished before final project 
completion may be realized. There must be only one final event repre- 
senting the job completion, and it must be possible to reach the final 
event from each node by following the arcs to successor nodes. 

In the usual PERT analysis as described in ref. 2 and 3, there are 
three estimates required for each activity time. They are used to 
describe a beta distribution and are a most likely amount of time 
labeled m, and optimistic*estimated on the time labeled a, and a 
pessimistic estimate labeled b. The mean time of completion (te) is 
computed by the formula: 


te = at+4mtb 
6 e 


The variance of the completion time associated with this mean is: 


v =( ba} 
6 


A mean and a variance is associated with each activity in the net- 
work. Then the network is analyzed “or the earliest expected project 
completion time and the variance associated with this estimate. 

A procedure and a computer program for solving a PERT network are 


given in ref. 2. The program was used as a model for PER-O-GRAPH 


eZ. 








described below. The procedure is explained using the network in 


Figure l. 





FIGURE 1. Example PERT Network 


The earliest expect time of occurrence, TE, ond the associated 
variance, V, for each event in the network is) computed by traversing 
the network from beginning events towards the final event. The activity 
preceding event 5 starts at time zero so the expected time of completion 
of event 5 is time 10 and the variance is 9. The activity connecting 
event 2-5 cannot begin until after the completion of event 5 or at time 
10. The earliest event 2 could occur on the basis of event 5 and its 
succeeding activity is 14 (10+4). However, before the TE for event 2 
can be determined the effect of other activities leading to event 2 must 
be considered. The earliest activity 2-4 could be expected to be 
completed is 27 (18+9). Therefore, activity 2-4 governs the TE of 
event 2, and the TE of event 2 is 27. The variance for event 2 is 
determined by adding the variance of its governing immediate predecessor 
event and the variance of the activity connecting that event to event 2. 


Therefore the variance associated with event 2 is 17 (1146). This 


process is continued until the TE and V for each event is determined. 


ie 








For the above example the TE for event 1, the final event, is 46. 
Therefore the earliest expected time of completion of the job is 46 
with a variance of 24. 

The next step in solving the PERT network is to determine the 
latest expected time at which an event can be allowed to occur without 
delaying the completion of the project. This value is labeled TL and 
is determined by traversing the network from right to left as follows. 
Event 2 can be allowed to occur as late as time 33 (46-13) and event 3 
can be allowed to occur as late as time 39 (46-7). Event 4 has two 
immediate successor events. The TL for event 4 is the smaller of the 
TL's determined by basing TL for event 4 on each of its in -ediate 
Successor events. Basing it first on event 2, the TL for event 4 is 
21 (33-9). Next basing it on event 3 the TL is 18 (39-21). Since this 
is the smaller of the quantities, the TL for event 4 is 18. The 
remainder of the network is then solved in this manner until the TL 
for each event is determined. 

Next the "Slack" of each event is determined. This is the difference 
between the latest time at which the event can be allowed to occur 
without delaying the overall project completion time and the earliest 
mime at which it car possibly occur. In formula form it is TL-TE. 
Events with a "Slack" equal to zero are critical events, and the path 
which joins the critical events is called the critical path. 

Figure 2 is the solution of the PERT network in Figure 1. The 
critica: path is the path along the arcs connecting events 4, 3, and l. 

The project completion times in the usual PERT analysis are assumed 


to approximate a normal distribution. In this example the project 
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completion times are assumed to be normal with mean 46 and standard 


AeviAtion Of 4.0. 


Event Numbers 


5 4 S 2 1 
Es 10 18 a7 27 46 
V 9 tt i ay 24 
TL ZS 18 ay he 46 
SLACK nie 0 0 6 0 


FIGURE 2. Solution for Example PERT Network 
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III. MONTE CARLO PERT ANALYSIS 


The use of Monte Carlo techniques offers many advantages over the 
usual PERT analysis described in the previous section. 

For a Monte Carlo analysis the network is drawn in the usual manner. 
Ho.;ever, instead of describing the distribution of each activity with 
the standard three parameters of m,a, and b, which assume a beta distri- 
bution, any distribution can be associated with an activity's duration 
time. This allows for greater flexibility and improved accuracy 
depending on how closely the distributions can be est” iated. 

A Monte Carlo analysis consists of a large number of simulations of 
the performance of the project. Each simulation has a possibility of 
estimating the time the actual project will take and the activities 
which will be critical. A simulation begins with the drawing of a 
random te from the distribution associated with each activity. The 
network is then solved in the usual manner except the variance is not 
computed. The project completion time and the identity of each arc 
which was on the critical path are saved. When the desired number of 
Simulations is completed, the overall project completion time distri- 
bution can be determined by plotting the individual completion times or 
by performing a mathematical analysis of the stored completion times. 
The mean and st-enrdard deviation of this distribution give a more 
accurate prediction of the project completion time than the usual PERT 
analysis since the distribution need not be normal as is assumed in 
the usual method. 

Rather Mian a prediction Of 42 Critical path. the eritrealtey ence « 


for each activity is given. This is simply the percentage of the 
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analyses for which each activity was on the critical path. This gives 
an indication of how much attention on activity should be given by 
management to ensure the earliest possible project completion. This 
can be very valuable, especially if the completion times on two or 


more paths are close. 
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IV. USES FOR MONTE CARLO PROGRAM 


The Monte Carlo PERT analysis program, PERT-O-GRAPH, is intended 
to demonstrate the way the critical path changes during a Monte Carlo 
analysis and to be an aid to a person learning this method to estimate 
project completion times. It can be used to solve small menagement 
problems, but due to the small buffer (4K) on the 2250 Display Unit at 
the Naval Postgraduate School for which this program was written, the 
arrays are limited so no more than fifty events and two hundred 
activities can be used. This number of events and activities could not 
be displayed simultaneously due to the buffer. Since most problems are 
larger than this, real life applications of the program are restricted. 

The displays are designed to demonstrate the effect of changing 
Mctivity duration cistributions upon the critic. lity indexes and the 
distribution of project completion times, and the effects of adding and 
deleting activities on the project completion times. 

PERT-O-GRAPH contains three main displays which are drawn on the 
2250 screen. The main display is that of the PERT network. The nodes 
and arcs are displayed wherever indicated by the user. At the bottom 
of this display are strings of text called light pen functions which 
enable the user to communicate with the program and indicate what 
function he wishes to perform next. If the user indicates he wishes a 
mean time analysis to bt performed, the mean time foi each activity is 
displayed on the arcs and the mean time critical path is indicated by a 
double line along the critical arcs. When performing a Monte Carlo 
analysis the user may indicate that he wishes the analysis to be per- 


formed in such a manner that after each single analysis the duration 
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time for each activity is displayed on its arc, and, when called, the 
eeitical path for that single analysis will be displayed. Jhis pro- 
cedure gives an instructor the option of having his students calculate 
the critical path by hand and then have the critical path displayed to 
check the answer. It also demonstrates how the critical path can vary 
depending upon the length of time required to complete the various 
activities. After running a complete analysis the indicated number of 
times this display has the criticality indexes drawn on each are in the 
network. This emphasizes the point that there is not necessarily only 
one critical path and that all activities must be supervised to varying 
degrees in order to realize the shortest project completion time. 

The second display is for entering the distribution of activity 
duration times. The user can indicate and enter either a cumulative or 
frequency distribution. Only «iscrete distributi: 1s can be entered 
with the choice of twenty values on both the X and Y axes. After 
entering the distribution and eae to the main ee the user 
may refer to any activity, and if a frequency distribution was entered, 
a plot of either the frequency of cumulative distribution associated 
with that activity will be displayed in the lower left of the main 
display. If a cumulative distribution was entered, only the cumulative 
distribution can be displayed. 

The third display is a plot of the distribution of the project 
completion times. The user designates the minimum time to be plotted, 
the width of the plot increments, and th: total number of increments. 

A plot of completion times takes up half of the display CRT. When more 
than one compli te analysis is run, the display shows th: plot of the 


most recent analysis and thuit of the previous analysis. This allows a 
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quick comparison of how altering the network changes the project 
completion times. 

In addition to the displays on the CRT there is a hard copy output 
for each completed analysis. This copy gives the sequential number of 
the analysis, how many iterations were made for that run, the mean and 
standard deviation of the project completion times, the upper limit of 
each increment, the per cent of total for each increment, and the 
cumulative percentage for each increment going from lowest to highest 


completion time recorded. 
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V. OPERATING INSTRUCTIONS FOR PERT-O-GRAPH 


The operating instructions. cover all the capabilities cf 7 ihiag: 
GRAPH in such a manner that the user need have no former experience 
with the IBM 2250 display unit and he need not have any computer pro- 
gramming ability. The use of the light pen, the 2250 keyboard, the 
programed function keys and the displays are explained. 

A PERT network entered into PERT-O-GRAPH has certain restrictions 
which must be met for successful execution. Node 1 must be the final 
event and it must be entered first. The rest of the noodes may be 
placed anywhere in the network, however, the internal execution of the 
program will be more efficient if the nodes are arran:.ed such that the 
lower numbered events are in the direction of the final evint. This 
arrangement is only an efficiency consideration and is not necessary 
for successful execuiion. A beginning event can have only one predeces- 
sor activity and the final,event cannot have any successors. All the 
events in the network cannot have more than four immediate predecessor 
events or four immediate successor events. However, a problem in which 
an event has more than four immediate predecessors or immediate succes- 
sors can be solved by using "dummy" events. The use of dummy events is 
illustrated in Figure 3. A maximum of fifty events can be entered and 
after event 50 appears no more cvents can be displayed. The maximum 


nunter of dif©erent distributions that may be defined is fifty. 
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FIGURE 3. Use of "dummy" Events 


fee LIGHT PEN FUNCTIONS 

Most of the user controlled operations are communicated to the 
computer by means of the light pen. The light pen operations which the 
user can call are displayed in Figure 4 on page 39 as they appear at 
the bottom of the main display. In order to call a light pen function 
the light pen must be placed over the desired text and then the light 
pen foot switch depressed. 

1. Enter Node 

Nodes are anteeea sequentially from 1. The first node entered 
will be 1, the second 2 etc. 

After using the light pen to detect the "Enter Node" function 
the user will be instructed to "Position light pen before pressing 
function key.'' The light pen should then be positioned over the place 
on the 2250 display screen where the user desires the next node to 
appear. Next the programed function key four lakeled "Locate L.P." 
should be pressed, and then the light pen foot switch depressed. When 
function key four is pressed the display will no longer be visible and 
the screen will be filled with W's. After depressing the light pen 
foot switch the PERT network will again be visible with the new node in 


the indicated position. All the nodes do not have to be entered before 


performing other functions. 


ee 





2. Delete Last Node 

The "Delete last node" light pen function causes the node just 
previously entered to be erased. If any function such as entering an 
arc has been performed since the last node was entered this option 
cannot be used, and the ''Delete node" function must be used to erase a 
node. If the "Delete last node" function is used the node deleted can 
be used again. This is not true if the "Delete node" function is used. 
The "Delete node" function is explained below. 

Souter Arc 

After indicating the desire to enter an arc, the user is instructed 
to 'Indicate node, then predecessor.'' By means of the light pen the 
user should indicate first the node, then the predecessor node between 
which the arc is to be drawn. After the arc appears the user has the 
@peion to delete the arc, ‘Delete last arc”, enter a new distribution 
for the new arc's activity time, "Enter distribution", and "Rad Dist 
from Card", or asse iate a ice ipeion which has already been entered 
with the new a “cane Uilsenibuetom as arc. 

If the "Delete last arc" option is selected with tie light pen, 
the arc will disappear and the original list of options will appear. 

After selecting the ‘Same distribution ac arc" option, the user 
indicates with the i 2ht pen which arc has the already defined distri- 
bution and the new are will share that di: tribution. 

The "Enter distribution" option causes the screen to be erased 
and the two op’ is "Cumulative distribution" and '°yvsquency distribu- 
tion" tc appear at the top center cf the display surface. The method 


for entering distributions is explained in B below. 
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4. Enter Beginning Arc 


The "Enter beginning arc' option must be used to enter the arc 
representing an activity with no predecessor event. After indicating 
this option the user is instructed to "Indicate beginning event." Then 
the three choices f.r entering distributions described under the Enter 
arc discussion are made available. 

5. Delete Node 

The "Delete node" light pen function allows the user to delete 
a node anywhere in the PERT network. After selecting this function the 
node to be deleted should then be indicated with the light pen, and it 
will be erased. The number on this node will not be reused and the next 
node entered will be sequenced as if no node were deleted. Since the 
final event must always be node 1, this option should not be used to 
delete node 1. After a node has been delcted any arcs leading to or 
from the deleted node must be deleted by using the "Delete arc" function. 

6. Delete Arc 

If the user wishes to delete an are for which a distribution of 
completion times has been entered, the ''Delete arc" function must be 
used. After choosing this function, the light pen is used to indicate 
the arc to be deleted, and that arc is erased from the network. 

7. Restart 

"Restart'' reinitializes the display. The PERT network displayed 
is completely erased as are all distributions of completion times pre- 
viously entered. The first node entered after a restart has ‘he number 
1. This fur-tion is usually used before a completely new network is 


entered for analysis. 
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8. Change Only Dist of Arc 
The "Change only dist of arc" function allows the user to 
change the distribution of an arc which is using a shared distribution 
without changing the shared distribution or to change the distribution 
of an arc which is not using a shared distribution. After selecting 
this function the user should indicate which arc is to have its distri- 
bution changed. Then the options to assign a new distribution or to 
share a distribution are presented. 
9. Change Common Dist of Arc 
The "Change common dist of arc'' function is used to change a 
distribution which is associated with more than one arc. After 
selecting this function, one of the arcs vwsing the common, or shared 
distribution should be indicated with the light pen. Then a new shared 
distribution may be entered, or a previously defined distribution can be 
associated with the arcs sharing the distribution by indicating which 
arc is associated with he dleelbaat distribution. This function may also 
be used to chare the distribution of an arc which is not sharing a 
distribution with other arcs, but ihe “Change only dist of arc" function 
is more efficient for this purpose. 
HOw Display Freq Dist of Are 
Ute) Display freq. disteon arCumtlUnctlommercmlced tomad sp lawsa 
frequency distribution associated with an arc. It can only be used if 
the distribution was entered as a frequency distribution. After 
choosing this function, the user indicates the are whose distribution 
is to be displayed, and the frequency distribution will appear in the 


lower left corner of the 2250 display. Pressing programed function key 3 


Zo 





labeled "Go on" erases the display of the distribution and allows 
further function calls. 
il. DispheysCumeDisiyoresrc 
The "Display Cum dist of arc" function is used to display the 
cumulative distribution associated with an arc. It does not matter 
whether the distribution was entered as a frequency or a cumulative 
distribution. After choosing this function, the are whose distribution 
is to be displayed should be indicated and the cumulative distribution 
will appear in the lower left corner of the 2250 display. To erase the 
display of the distribution and allow further function calis the pro- 
gramed function key 3 labeled "Go on" must be pressed. 
ee Out 
The "Quit" light pen function ends the program and disables the 
2250 display unit. 
13. Execute 
The "Execute" light Pemeeinc cies indicates that the PERT network 
is complete and ready to be analysed. After selecting this light pen 
function the user has no means of changing the network until the 


analysis is completed. 


B. PROGRAMMED FUNCTION KEYS 

The IBM 2250 Display Unit has associated with it a programed function 
keyboard. There are 32 keys numbered from 0 to 31. In PERT-O-GRAPH, 
keys 0 through 6 are used. The lights inside these seven keys are on 
only when they may be used, however the user must be careful to press 
the correct key depending on what part of the program is currently being 


executed. Pressing a key causes an interrupt and a message to be 
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communicated to the program, and pressing the wrong key may cause 
unpredictable results. 
1.9 End Dist.(P aie heya 
This P.F. key is used when entering a distribution of activity 
completion times. It indicates that the distribution has been 
completely entered. Before pressing this key the user should check that 
the distribution is correct, since after it is depressed the distribu- 
tion is entered into the program and one of the functions to change the 
distribution of an are would have to be called to correct the 
distribution. 
ga) End cf Char (P.F. Key 1) 
pcm ea eey el Suusedeat Let cnberun: sedarachere: Dyotnewalpna. 
numeric keyboard. Pressing this key will ee cursor from the 


character string and enter the character string into the program. 


Bee Restart Plot (P.F. Key 2) 


a 


This P.F. key can be used when entering the plot of a distribu- 
tion of activity completion times. Jt causes the labels on the X axis 
and the plot to be erased. It will be recognized only after the first 
point is plotted. 

4. Go On (P.F. Key 3) 

This P.F. key comaunicates to PERT-O-GRAPH that the user wishes 
the program to execute the next step in the program. The result depends 
upon the step presently in progress. 

a. Singl Step Execution 

When pe: forming a Mo. te Carlo analysis in single step mode, 


the te for each are is displayed when the single analysis is completed. 
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Pressing this key after the te's are displayed will cause the critical 
path to be drawn. 

b. Criticality Indexes Displayed 

Pressing P.F. key 3 when the criticality indexes are 
displayed after a complete Monte Carlo analysis will cause the PERT 
network to disappear and the plot of the project completion times to be 
drawn. 

c. Distribution of Project Completion Times Displayed 

Pressing P.F. key 3 when the distribution of project 
completion times is display’d will cause the distribution to be erased 
from the screen, the PERT network to appear, and will allow alteration 
of the networ: . 

d. Mean Tim. Calculation 

After a mean time calculation, the te's and critical path 
are displayed. Pressing P.F. key 3 at this time will cause the te's 
mmaecritical path to be eee allow the user to enter the number of 
analyses to be pericrmed in the Monte Carlo analysis. 

me Locate |..P. (P.F. Key 4) 

After indicating that a new node is to be entered, pressing P.F. 
key 4 will cause the 2250 display area to be filled with W's. This is 
necessary for the light pen to detect the coordinates where the new node 
is to bc drawn. 

Gee Step (P.F. Key 5) 
This key is used when performing a Monte Carlo aualysis in the 


Single execution mode. After the critical path is displayed pressing 


this key will cause the newt analysis to be performed. 
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/ Conti nuewm lomneces (P.F. Key 6) 
After a complete Monte Carlo analysis has been performed and 
the plot of project completion times is displayed, this key is used to 
indicate that the same network is to be solved an additional number of 


times. 


C. ALPHANUMERIC KEYBOARD 

The alphanumeric keyboard attached to the 2250 Display Unit allows 
the user to enter characters into the display. In order for a character 
to be entered into the display a cursor must be in the display. The 
cursor is inserted autoatically by PERT-O-GRAPH een information is to 
be entered via this keyboard. After a character is ente :d, the cursor 
moves to the next position where a character can be entered. 

The ADVANCE key moves the cursor one position without changing the 
character above the cursor. Pressing P.F. key 1] will remove the cursor 
and enter the text into PERT-O-GRAPH. Hitting the ALT and CANCEL keys 
Simultaneously will cause the display to disappear and the text "*TER- 
MINATE", "*DUMP", and '*RESUME" to appear. Detecting on "*TERMINATE" 
and "*DUMP" with the light pen will cause the program to terminate, 


and "*RESUME" will return the display and execution will continue. 


D. DRAWING A DISTRIBUTION 

After entering an arc or indicating the desire to change the distri- 
bution of an arc and then indicating a new distribution is to be entered 
by the user the 2250 screen will be erased and the text "Cumu‘ative 
distribution" and "Frequency distribution" will appear at the top center 
ee the display. The lig't pen is used to indicate the type of distri— 


bution to be entered. This will cause an X and Y axis to appear. 
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Each axis is divided into twenty segments by strokes along the axis. 
mie Y axis is ljabeled from0.0 to.1.0- 
io Oo bapeling Pie as 

The extremes on the X axis are entered by the user. A cursor 
first appears to the left of the axis where the left extreme is 
entered. The extremes are entered by using the 2250 display unit key- 
board. After entering the left extreme, pressing programed function 
key 1 labeled "End of char." will cause the cursor to appear to the 
right of the X axis. The right extreme is entered in the same manner. 
The axis will bce labeled along its length with increments of one-tenth 
of the difference between the extremes. It is not necess: -y to enter 
.he most optimistic and most pessimistic times zs the exti: -es. The 
extremes should be chosen o the axis will be labeled with points that 
are to be used in the plot. As an example, if the user wished to plot 
points between 7 and 13 he might enter the extremes of 5 and 15. This 
would cause 5.0, 5.5, 6.0, 6.5, etc. to be entered along the axis. If 
7 and 13 were entered as the extremes, the labels would be 7.0, 7.6, 8.2, 
9.4, etc. The inaximum value that can be entered as either extreme on 
the X axis is 9999.9. Ali values entered must contain a decimal point 
and may contain one digit to the right of the point. 

2. Drawing a Cumulative Distribution 

To enter the starting point of a distribution the light pen is 
used to indicate the vertical hash mark on the X axis representing the 
most optimistic estimate. The distribution is drawn by first indicating 
the X value of the point by detecting one of the hash marks on the X 
axis and then indicating the Y value by detecting one of the hash marks 


on the Y axis. PERT-O-GRA!:! will draw lines between succeeding points. 


30 








Any discrete cumulative distribution of twenty points or less can be 
drawn. The last point entered must have a Y value of 1.0. When 
finished entering the distribution, programed function key 0 labeled 
"End dist'’ should be pressed. The distribution will then be entered, 
and the display of the PERT network will reappear. 
3. Drawing a Frequency Distribution 

To enter the starting point of a distribution the light pen is 
used to indicate the vertical hash mark on the X axis representing the 
most optimistic estimate. To draw he distribution first indicate the 
X value of the point by detecting one of the hash marks on the X axis 
and then indicate the position of the Y value by detecting one of the 
hash marks on the Y axis. The values on the Y axis can be thought of 
as positions instead of values since PERT-O-GR/‘PH will draw lines 
between succeeding points and internally adjust the values so the area 
under the finished plot is equal to 1.0. When finished entering the 
distribution, programed function key 0 labeled "End dist" should be 
pressed. The distribution will then be entered, and the display of the 
PERT network will reappear. 

eeEsor Correction 

If an error is made when entering a distribution, the plot and 
the X axis labels can be erased by pressing programed function key 2 
labeled "Restart plot''. This key can be used any time after the first 


point is entered. 


E. ENTERING A DISTRIBUTION FROM PUNCHED CARDS 
A distribution may also be entered from punched cards by detecting 


mie inter dist from card’ light pen function. If distributions are to 
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be entered using this option, the punched cards must be entered at the 
end of the program deck prior to submitting the program for execution. 
Only cumulative distributions can be entered through this function. 

The first card for each distribution must contain the number of 
points used to describe the distribution. This number must be an 
integer and it must be placed in the first six columns of a card. The 
last digit must be in column six as any trailing blanks will be interpre- 
ted as zeroes. 

The points in the distribution are entered on separate cards. The 
first number on a card is read as the X value. It must be entered in 
the first six columns of a card. If no decimal point is entered the 
defauli format of F6.1 will be assumed. The second number on a card is 
read as the Y value of the point and it must be ent.red in card columns 
seven through twelve. If no decimal point is entered, the default 
format of F 6.3 will be assumed. 

The user should be careful to place the distributions in the order 


in which they are to be entered as they cannot be read out of order. 


Once a distribution is read from cards it cannot be read again. 


Pee EXECULION 

After indicating that an execution is desired the original list of 
options is erased from the screen and the two options ''Mean time 
calculation", and "Execute-Analyses" appear. The mean time calculation 
will perform only one analysis and display the mean -ime of project 
compli tion using thc mea: “ime for each activity. This is similar to 
the usual PERT analysis. !.y indicating the number of analyses to be 


made the user indicates that a Monte Carlo analysis is to be made. 
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1. /Mean time, Catlewlacion 

When a mean time calculation is selected with the light pen the 
mean time is written to the right of "Mean time calculation" and the 
mean time of each activity's completion time is displayed on its arc. 
Also the mean time critical path is indicated by a double line along 
the arcs on the critical path in the PERT network. The user may now 
indicate the number of analyses to be performed in the Monte Carlo 
analysis. 

2. Monte Carlo Analysis 

To begin a Monte Carlo analysis enter the number cf analyses to 
be run in "Execute —- Analyses" with the 2250 keyboard. The largest 
number that can be entered is 9999. If less tha: 1000 analyses are 
desired, this number must end with a decimal point. After this number 
is entered, function key 1 should be pressed, and then the user will be 
direc: . to enter "Minimum time." This is the minimum time that will 
be plott:d in the etriomtion of project completion times. The largest 
number that can be entered here is 999999.9 and a decimal point must be 
included. Next the user is directed to enter the "Width of increments.” 
This is the distance between points in the plot of project completion 
times. This number also can be no larger than 999999.9 and must contain 
a decimal peint. After again pressing function key 1, the user is 


@epeeted to enter "No. of increments." 


This is the number of points to 
be plotted in the distributi. of project completion times. The 
maximum nur er of points i 50, and the numb: r must end with a decimal 
poant, 

After entering this information the user must make the choice of 


"Single execution" or "Complete analysis." 
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3. Slag Lee ceut Tom 

If the "Single execution" mode is chosen with the light pen, 
one analysis will be performed, and the te picked for each activity 
fil) be written on that activity'’s arc. Next, prosran@@) function key 3 
should be pressed, and the critical path will be displayed by a double 
line along the arcs of the critical path. If another single execution 
is desired, it may be performed by pressing function key 5 labeled 
"Step" or by detecting "Single execution" with the light pen. After 
any Single execution the user may go to the complete analysis mode by 
indicating "Complete analysis" with the light pen. | 

4, Complete Analysis 

When the "Complete analysis'' mode is selected, the nunxber of 
Monte Carlo analyses indicated are perforred, and when they are completed 
the criticality index fer each are is written along each arc. This 
ind“ cates the percentage of the time an arc was on the critical path. 
After the criticality indexes are displayed, pressing programed function 
key 3 will cause the PERT network to be erased and a plot of the project 
completion times to be displayed. The current plot fills half of the 
Screen and the previous plot fills the other half. 

At this point, the user has two options. He can press function 
key 3 ("Go On") and the prograr. will return to the point where modifica- 
tion of the PERT network is allowed. If function key 6 labeled "Continue, 
no reset" is pressed the user will be directed to enter the number of 
additional annlyses to be performec end then to indicate step execution 
or complete execution. The "Continue, no reset" function does not reset 
any of the data already compiled and allows for additional executions to 


em 20th the curve of project completion times. 
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G. SUBMITTING A RUN 

PERT-O-GRAPH is complete by itself and can be run by a user with 
no additions except a JOB card unless he wishes to read distributions 
from cards. The green JOB card is placed at the beginning of the 
program deck and it must contain the phrase TYPERUN=HOLD. Cards 
containing distributions of activity duration times must be placed 
immediately before the orange # card which is at the end of the 


program deck. 
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VI. SAMPI PROBLEM 


The problem described in this section is to determine the distri- 
butions of activity duration times necessary to cause the plot of 
project completion times for the network shown in Figure 4 on page 39 
to have a mean below 50.0. The network represcnts the jobs in the 
project which must be completed and the order in which the jobs must 
be started. Node 1 was entered first, then the remaining nodes and 
arcs were entered. 

Some of the distributions of activity duration times were entered 
from punched cards and some were entered by drawing the distributions 
on the display unit. The distribution for arc 5-8 was drawn, and is 
shown in : igure 5 on page 40. Before the begii. ing arc entering node 
7 was enterec, the distribution for are 8-5 was displayed as shown in 
Figure 6 on page 41. Since this was the sams as the distribution 
desired for the beginning arc entering node 7 this arc was designated 
to share the distribution of arc 5-8. 

There were six distributions entered for the twelve arcs. The 
points in the distributions and the arcs which were associated with 
these distributions are given in Table 1 on meee SON 

After the complete network was entered, a mean time calculation was 
executed. This resulted in a mean time of 51.4 and a critical path 
along the arcs connecting nodes 8,6,3, and l. 

Next a Monte Carlo analysis of 1,000 cxecutions was spe’ified. 
Using the mean time calculation as a basis, a minimum time of 30.0, a 
wi' “h of increment of 1.5, and the number of increments 30 were entered 


a 


for the plot of projcct completion times. When this analysis was 
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completed, the mean time of the plot of the project completion times 
was 54.9 with a standard deviation of 4.0. It was then run again for 
another 1000 executions. 

The criticality indexes bas Jd on these 2000 executions are shown 
in Figure 7 on page 42 . The mean time of project completion times 
resultine from those 2000 analyses was 54.9 with a standard deviation 
of 4.1. Since these values were so close to the values obtained from 
1000 analyses, additional executions to smo: th the curve of completion 
times and to calculate the criticality indexes were felt unnecessary. 

The estimated mean tiie of project completion using a Monte Carol 
analysis was 54.9 as compared to a mean time of 51./ using a mean time 
calculation which is similar to the usual PERT analysis. This is the 
Same type o” results as were described by Van Slyke in ref. 3. 

Next an attempt was made to lower the mean time of the project 
completion times. The arcs with the highest criticality inde:es 
seemed to be the logical ones to alter. Arc 1 - 3 had a criticality 
index of ./3 so this arc was changed first. A uniform distribution 
between the times of 5.0 and 12.0 was entered for this arc's distribu- 
tion of duration times. 2000 analyses were specified and the same 
parameters for the plot of project completion times were entered. This 
is shown in Figure 8 on page 43. 

This analysis was first run in the single step mode so the changing 
of the critical pati could be observed. The resulting critical path 
from two individual an. lyses re shown in Figuics 9 and 10 on pa; es 44 
and 45, 


The analysis was then run in the complete analysis mode to the 


finish. This resulted in a plot of project completion time with a mean 
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of 53.5 and a standard deviation of 4.2. The plots of the two complete 
analyses are shown in Figure 11 with the old plot at the top of the 
display and the new plot at the bottom. 

Another attempt was then made to lower the mean of the distribution 
of the project completion times. The arc with the highest criticality 
index in the last analysis was arc 4 - 7 with an index of .65. The 
distribution of duration times for this arc was changed to a cumulative 
distribution uniform over the times from 7 to 15. The network was then 
analysed 2000 times with a resulting distribution of project completion 
times with a mean of 49.1 and a standard deviation of 4.0. 

This problem illustrates how the user can experiment with a PERT 
network using PERT-O-GRAPH. He can see how the critical path changes 
with differen. activity durations, and he can observe the effect of 


altering the network on the distribution of project completion times. 
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PELT Network Analyzed in Example Problem 
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FIGURE 5. Cumulative Distribution drawn by light pen 
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*IGURE 6. Display of Cumulative Distribution associated 
with ar Arc. 
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FIGURE 7. Display of Criticality Indexes 
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FIGURE 8. Display of Param.ters for Plot of Distribution 
of Project Completion Times. 
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. FIGURE 9. Display of Network with critical path shown 
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Display of Network with Critical Path changed 


FIGURE 10. 
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- FIGURE 11. Display of Plots of Distributions of Project 
Completion times 
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PERT-O-GRAPI] is written for an IBM 360/67 with an IBM Display Unit 
Model 1. The Display Unit 1s deseribed in detail an ref. 4. The 
proyram is written in FORTRAN IV, described in ref. 7 tO: Gun on a 
machine with 0S/360 release 18 and an IBM Graphics Subroutine Package 
described in ref. 6. 

The program for the PERT analysis is modeled after a PERT program 


given in ref. 1. The model program has been expanded to handle nodes 


with four immediate successor and four immediate predecessor events. 


A. MAIN PROGRAM 
The main program controls the 2250 displays, the calling of sub- 
routines, the filling of the various arrays, deteriining the critical 
path and completion time of the PERT network, and the printing of the 
hard copy output. 
ee the 2250 Displays 
The display fox the PERT network is entered into its own graphic 
data set. This data set is equivalenced with two other data sets. One 
of these data sets is used when entering a distribution and one is used 
in the ple of project completion times. When one of these other data 
sets is called, the data set with the network is overlayed in the 2250 
buffer and a copy is retained in the 360's core. When the data set 
containing the network is recalled, it overlays an equivalent data set 
in the buifer and reappears. 
The nodes in the network have a correlation value assigned so 
they can be identified by light pen detections. The correlation value 


assigncd is the same as the number of the node. The arcs in the 
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network are correlated with a number which identifies the two nodes 
between which the arc is drawn. For exam,le, an arc correlation value 
of 0310 identifies the arc going from node 3 to predecessor node 10. 

The text which appears on the screen for light pen functions 
and instructions to the user is placed in its own data set. This data 
set is completely erased when the program jumps to another section of 
code. The reason this text is not written just once and then placed in 
include or omit status when needed is that it was found that the text 
took up too much buffer space when used in this manner. By completely 
mesetting the dita set, the buffer space used Ly this data set is 
freed for other uses. 

The display for entering a distribution is composed of two 
graphic data sets. One data set contains the axes, the labels on the. 
axes, and the plot. The oth. r contains only the hash marks on the X 
and Y axes. The light pen is enabled to detect only the data set which 
contains the hash marks. Only ai characters can be entered for the 
left and right extremes. This allows a space equal to the size of one 
half character between strings of numbers after the entire X axis is 
labeled. When the user enters the X coordinate of a point, the program 
saves the absolute X coordinate for plotting the point. The time value 
associated with this point is computed by scaling. The coordinate and 
time value for the Y axis are handled similarly. The same display is 
used to enter a frequency or a cumulative distribution. 

The piot of the distribution of project completion times 
involves the use of two data sets. After two complete analyses, both 
data sets contain a plot. After the third analysis, the oldest plot is 


erased and the new plot is entered in the cleared data set. In this 
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manner the current and last distribution are displa;cd simultaneously. 


7 


The extremes on the X axis are determined by the valucs the user cnters 


tina 


for "Minimum Time", "Width of Increment," 


and "Number of Increments." 
The smallest increment includes the number of completions which fall 
in that increment as well as those that are less than the minimum value. 
The largest valued increment includes the number of completions which 
fall in that increment as well as those that are greater tha: the right 
extreme. The Y axis is labeled from 0.0 to the maximum percentage of 
completions any increment contains. The plot is a series of straight 
lines. The X coordinates for the end of each line plotted are the 
centers of the increments except for the two extremes. ‘lhe plot starts 
at the Y axis < J ends with the right extreme. The value of the Y 
coordinate represents the percentage of comple! ion times in each 
increment. 
Baeearrays 

There are three main arrays in the program. The EVT array, 
dimensioned 50 by 20, contains the data necessary for the PERT analysis. 
The IPT array, dimensioned 50 by 4, is linked to the EVT array and 
contains pointers to the array containing the distributions of activity 
duration times. The DIS array is dimensioned 50 by 20 by 3 and contains 
the distributions of activity duration times. 

The first subscript of the EVT array indicates the number of 
the node. The data contained in the EVT array is explained in Table 1 
on pag. 20. 

The first subscript of the IPT array indicates the number of 


the node. The second subscript indicates the location of a pointer in 


the array for a particular arc's distribution of duration times. 
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EVT(N,1): 


EVT(N,2): 
EVT(N,3>: 
EVE (N. 4): 
EVT(N,5): 
evil (N.6) * 
VEC, 7): 
EVT(N, 8): 


BVI(N,9): 


mvi(N, 10) : 


Ba (N, lL): 
Mvr(N, 12): 
myrCN, 13): 
EVI(N,14): 
EVT(N,15): 
EVvE(N, 16): 
EVT(N,1.7): 
Hy@CN. 18): 
EVT(N,19): 


BVI(N, 20) ; 


Number 


of event N's first I. P. (immediate 


Predecessor) 


The te 


Number 


The te 


Number 


The te 


Number 


The te 


associated with N's first I. P. 
of event N's second I. P. 
associated with N's second i. P. 
Of GYent N Sethira 1. e. 
associated with N's thasrd =. 2: 
Of event ce sourth fol. 


aSSoOciatea with oN S fourth 2». e. 


TE of event N 


Number 


of event N's first I. S. (Immediate 


Successor) 


The te 


Number 


The te 


Number 


The te 


Number 


The te 


AS 6 OG eclawW isehie Ni oeuk tic tu Ge 
of eve it N's second I. S. 
associated with N's second I. S. 
of event N's third I. S. 
associated with N's third I. S. 
of event N's fourth I. S. 


asSociatcduwitheN s fourth, £. os. 


TL of event N 


SLACK of event N 


O if event current, 1 if event deleted 


DPABIGE 1. 


Definition of EVT Array 
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For example IPT(N,1) contains a pointer to the distribution of duration 
times for event N's first immediate predecessor, IPT(N,2) points to the 
distribution for N's second immediate predecessor, etc. 

The DIS array contains the distributions of activity duration 
times. All the distributions are stored as cumulative distributions 
even if they are entered as frequency distributions. The first sub- 
script is the number of the distribution and is pointed to by the IPT 
array. The second subscript is the number of the point in the distri- 
bution. A distribution must be described with twenty points or fewer. 
DIS (1,J,1) contains the cumulative distribution's Y value for the 
Jth point of the Ith distribution; DIS (1,J,2) contains the corresponding 
X value; and DIS (1,J,3) contains the corresponding Y value for the 
frequency distribution if thc user enters the distribution as a 
{ 2quency distribution. If the user enters a cumulative distribution, 
the cells for the Y value for the frequency distribution are set to 
ZELO . 

When an arc is entered, the identity of the predecessor is 
entered into the EVT array. Then the pointer to the distribution of 
activity duration times is entered into the IPT array. 

The ITCRT and ICRIT arrays are dimensioned 50 by 4. The ITCRT 
array kee:.s the information on which arcs are on the critical path when 
doing a mean time calculation or a Monte Carlo analysis in the single 
step mode. The ICKIT array kceps a count on how many times each arc 
was on the critical path for the criticality ind.« calculation. The 
IHIS array is a vector of length 50 which contains the total number of 


completions which fall into each increment. The LOC array is dimensioned 


50 by 2 and contains the X and Y coordinates of each node. The PLOT 
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array is a 50 by 2 array of points to be plotted by the DISPLY sub- 
mOouLine..— The P arreyeis a vector of length 150 used as a stack of 
pointers in the analysis of the PERT network. 

3. Analyzing the PERT Network 

If the user desires a 'Mean Time Calculation" to be performed, 
he must indicate this before entering the number of iterations to be 
performed in a Monte Carlo analy’is. 

Before determining the mean time critical path and completion 
time, the mean duration of each arc is calculated. This is determined 
by calling subroutine MEAN for each arc and then entering the mean time 
for Gach arc's distributien into the EVI array as that arc'’s te. The 
analysis is then run once and the mean time for each arc's duration, 
the mean time for completion, and the mean time critical path are 
displayed. ‘the user may then chooose a Monte Carlo analysis. 

In a Monte Carlo analysis, before determining the eriered 
path and completion time of a network, a te is determined for each arc. 
This is accomplished by generating a random number between 0.0 and 1.0 
in the RAND subroutine, entering the SETTE subroutine with this number, 
and calculating a te from the associated ‘istribution for an arc. The 
te is then entered in the ap»ropriate place in the EVT array. This 
procedure is repeat«? for each arc in the network. After all the te's 
are entered, the TE for each node is computed. The program then 
calculates the TL for each node. Next the "SLACK" for each node is 
calculated by su tracting the TE from the TL for each node. 

ine prog m next determines which arcs are on the critical 
path. This is done by first checking each node's slack. If the slack 


is zero, it checks each are that entered the node, and if the TL of the 
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node minus the te associated with the arc equals the TL of the 
predecessor, and the predecessor also has a slac’. of zero, then that 
mec 1S on the critical path. The critical arcs are tabullated angie 
ICRIT array for a complete analysis and in the ITCRT array for display 
of the single step critical path. This process is repeated the number 


of times specified by the user. 


ieee oUBROULINES 

The function of each of the subroutines in PERT-O-GRAPH and 
descriptions oi their e:. cution are included in this section. 

ie ouproutine DIST 

Thiseeueroutime epters a distribution of an activity s duraciom 

times into the D.S array. If a cumulative distribution is drawn on 
the 2250, it is entered into the array exactly as plotted. Ifa 
frequency distribution is drawn, the main program changes it into a 
cumulative distribution. The user is cautioned that the cumulative 
distribution resulting from this conversion is represented internally 
as a series of straight lines, the ae number of straight lines as was 
used to draw the frequency distribution. If a greater number of lines 
is used to draw the frequency distribution, the resulting cumulative 
error will be reduced. For example, a triangular distribution drawn 
with two lines becomes a straight line cumulative distribution. However, 
if may short li segments are used to draw the same frequency 
distribution, the resulting cumulative will more nearly resemble the 
correct curved line cumulative distribution. Figure 12 illustrates 


tais result. 
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An unavoidable error is caused by the 2250 light pen detections. 
The 2250 rounds tne raster unit coordinates of a detection so that 
when the coordinates are converted to values on the X and Y axes, 
slight errors may result. 
2. Subr itine MEAN 
This subroutine calculates the mean of the distribution of an 
activity's duration times. The mean is determincd by calculating the 
average X value under each line segment of the cumulative distribution, 
weighing these values according to the amount of the Y axis spanned by 
each line seginent and taking the sum of these weighted values. This 
sum is returned as the mean of the distributio:. 
3. Subreutine RAND 
This subroutine is a random number genereior provided in ref. 7. 
It starts with a seed of 65549 and generates random numbers between the 
wemumes of 0.0 and 1.0. 
fee subroutine SEITE 
This subroutine enters a distribution of an activity's duration 
times with a random number and determines a te for the activity. The Y 
value of the cumulative distribution is searched sequentially until 
value: brecketing the random number are found. Then the X value cor- 
responding to these Y valucs is interpolated and the resulting X value 
is returned as the te. This is the inverse transiormation method of 
generating stochastic variates as described in ref. 8. 
we SuenrOuLdne EL iias 
This subroutine fills the PLOT array for plotting t!e distribu- 
tion of project completion times. It labels the plot of completion 


times, determines where the plot is to be placed in the display, and 


calls the DISPLY subrout. ne. 
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FIGURE 12. Drawing a Frequency Distribution 
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b.. Subsoueine Dre iy 
This subroutine is called to display the plot of project comple- 
tion times and to display the distribution of an arc which has already 
been entered. It determines how many labels will fit on the X and Y 
axis, labels the axes, and plots the PLOT array. 
i oubroutine CRIT 
This subroutine uses the ITCRT array and accentuates the arcs 
on the critical path. It is called after a mean time calculation and 
when performing a Monte Carlo analysis in the single step mode. Two 
extra lines are drawn along the arcs on the critical path. One line is 
displac..:i four raster units to the left and four up from the arc. The 
second line is displaced foux raster units to the right and four up 
meom the arc. 
8. Subroutine MARK 
This subroutine generates text above the arcs in the network. 
It displays the mean time te when doing a mean time calculation, the 


te when doing a Monte Carlo analyses on the single step mode, and the 


criticality indexes after the completion of a Monte Carlo analysis. 


a0 





VIII. CONCLUDING REMARKS 


The program described in this report allows a user to draw a PERT 
network on an IBM 2250 graphic display unit. The user can also draw 
distributions of activity duration times or he can enter them from 
punched cards. The network is then analyzed by the computer and the 
results of the analysis are displayed on the graphic display unit. 

One of the uses of the program is to demonstrate for typical 
problems some of the differences between the usual PERT analysis and a 
Monte Carlo PERT analysis. In the usua? PERT analysis the network is 
solved only once and the result is one critical path and a project 
completion time based on the mean times of activity durations. This 
is what is done when a network is analyzed for a 'Mean Time Calculation" 
in tis program. In a Monte Cz lo analysis the project is simulated a 
large number of times to produce a distribution of project completion 
times. For each simulation the activity duration times may be 
different, with resulting differences in the critical path. A 
criticality index is associated with each activity as a measure of its 
chances of being on the critical path. A Monte Carlo analysis usually 
estimates a project completion time distribution with a mean time 
greater than that of the usual PERT analysis. Also, it usually shows 
that there are many potential critical paths. 

Another use of the program is to allow the user to experiment with 
making changes in the network. Nodes and arcs can be added or deleted 
and the distributions of activity durations can be changed. The user 
can then analyze the changed network and see what effect the changes 
have on the criticality indexes and on the distribution of project 
completion times. 
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An improvement that could be made to the program would be to modify 
it to accept an entire network, with associated distributions, from 
punched cards. If the same network is to be analyzed and experimented 
on many times this would save the user the time required to reenter 
the network by light pen before each experiment. 

Another improvement which could be made would be to add. a sub- 
routine which whenever called would plot whatever is displayed on the 
graphics unit on a Calcomp Plotter. This would give the user the 
capability of getting a hard copy output of the network, the plots of 


the distributions or any other display which he might find useful. 
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